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Abstract 

The rapid rise of global temperature that began about 1975 continues at a mean rate of about 
0.18°C/decade, with the current annual temperature exceeding +1.25°C relative to 1880-1920. 
Global temperature has just reached a level similar to the mean level in the prior interglacial 
(Eemian) period, when sea level was several meters higher than today, and, if it long remains at 
this level, slow amplifying feedbacks will lead to greater climate change and consequences. The 
growth rate of climate forcing due to human-caused greenhouse gases (GHGs) increased over 
20% in the past decade mainly due to resurging growth of atmospheric CH 4 , thus making it 
increasingly difficult to achieve targets such as limiting global warming to 1.5°C or reducing 
atmospheric C0 2 below 350 ppm. Such targets now require “negative emissions”, i.e., extraction 
of C0 2 from the atmosphere. If rapid phasedown of fossil fuel emissions begins soon, most of 
the necessary C0 2 extraction can take place via improved agricultural and forestry practices, 
including reforestation and steps to improve soil fertility and increase its carbon content. In this 
case, the magnitude and duration of global temperature excursion above the natural range of the 
current interglacial (Holocene) could be limited and irreversible climate impacts could be 
minimized. In contrast, continued high fossil fuel emissions by the current generation would 
place a burden on young people to undertake massive technological C0 2 extraction, if they are to 
limit climate change. Proposed methods of extraction such as bioenergy with carbon capture and 
storage (BECCS) or air capture of C0 2 imply minimal estimated costs of 104-570 trillion dollars 
this century, with large risks and uncertain feasibility. Continued high fossil fuel emissions 
unarguably sentences young people to either a massive, possibly implausible cleanup or growing 
deleterious climate impacts or both, scenarios that should provide both incentive and obligation 
for governments to alter energy policies without further delay. 
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Fig. 1 . Fossil fuel (and cement manufacture) C0 2 emissions based on Boden et al (2016) with BP data 
used to infer 2014-2015 estimates. Europe/Eurasia is Turkey plus Boden et al categories Western Europe 
and Centrally Planned Europe. Asia Pacific is sum of Centrally Planned Asia, Far East and Oceania. 
Middle East is Boden et al Middle East less Turkey. Russia is Russian Federation since 1992 and 0.6 of 
USSR in 1850-1991. Ships/Air is sum of bunker fuels of all nations. Can+Aus is Canada+ Australia. 


1. Introduction 

Almost all nations agree on the need to limit fossil fuel emissions to avoid dangerous human- 
made climate change, as formalized in the 1992 Framework Convention on Climate Change 
(UNFCCC 1992). The Paris Agreement (2015) seeks to limit global warming to well below 2°C 
relative to pre-industrial levels, with an aspirational goal of staying below 1.5°C.We advocate a 
stricter goal, based on restoring Earth’s energy balance and limiting the period when global 
temperature is above the range of the Holocene; temperature stability of the Holocene has 
allowed sea level to be stable for the past several millennia in which civilization developed. This 
goal leads to a C0 2 target of 350 ppm, which can be adjusted as C0 2 declines and empirical 
climate data accumulates (Hansen et al 2008, 2013, 2016). Either target, 1.5°C or 350 ppm, 
requires rapid phasedown of fossil fuel emissions. 

Despite widespread recognition of the risks posed by climate change, global fossil fuel 
emissions continue at a high rate that tends to make these targets increasingly improbable. 
Emissions are growing rapidly in emerging economies; while growth slowed in China in the past 
two years, emissions remain high (Fig. 1). The Kyoto Protocol (1997), a policy instrument of the 
Framework Convention (UNFCCC 1992), spurred emission reductions in some nations, and the 
collapse of the Soviet Union caused a large decrease of emissions by Russia (Fig. lb). However, 
growth of international ship and air emissions (Fig. lb) largely offset these reductions and the 
growth rate of global emissions actually accelerated from 1.5%/year in 1973-2000 to ~2. 5%/year 
after 2000 (Fig. Al). China is now the largest source of fossil fuel emissions, followed by the 
U.S. and India, but on a per capita historical basis the U.S. is 10 times more accountable than 
China and 25 times more accountable than India for the increase of atmospheric C0 2 above its 
preindustrial level (Hansen and Sato 2016). Tabular data for Figs. 1 and Al are available on the 
web page www.columbia.edu/~mhsl 19/Burden. 

In response to this situation, a lawsuit [Juliana et al vs United States 2016, hereafter J et al 
vs US 2016] was filed against the United States asking the U.S. District Court, District of 
Oregon, to require the U.S. government to produce a plan to rapidly reduce emissions. The suit 
requests that the plan reduce emissions at the 6%/year rate that Hansen et al (2013) estimated as 
the requirement for lowering atmospheric C0 2 to a level of 350 ppm. At a hearing in Eugene 
Oregon on 9 March 2016 the United States and three interveners (American Petroleum Institute, 
National Association of Manufacturers, and the American Fuels and Petrochemical Association) 
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asked the Court to dismiss the case, in part based on the argument that the requested rate of fossil 
fuel emissions reduction was implausible. Magistrate Judge Coffin stated that he was “troubled” 
85 by the severity of the requested emissions reduction rate, but he also noted that some of the 

alleged climate change consequences, if accurate, could be considered “beyond the pale”, and he 
rejected the motion to dismiss the case. Judge Coffin’s ruling must be certified by a second 
judge, after which the case can proceed to trial. It is anticipated that the plausibility of achieving 
the emission reductions needed to stabilize climate will be a central issue at the trial. 

90 Urgency of initiating emissions reductions is well recognized (IPCC 2013, 2014; 

Huntingford et al 2012; Friedlingstein et al 2014; Rogelj et al 2016a) and was stressed in the 
paper that the lawsuit J et al vs US (2016) uses to prescribe an emissions reduction scenario 
(Hansen et al 2013). The climate research community also realizes that the goal to keep global 
warming less than 1.5°C probably requires negative net CO 2 emissions later this century if high 
95 global emissions continue in the near-term (Fuss et al 2014; Anderson 2015; Rogelj et al 2016b; 
Sanderson et al 2016). The Intergovernmental Panel on Climate Change (IPCC) reports (IPCC 
2013, 2014) do not address environmental and ecological feasibility and impacts of large-scale 
CO 2 removal, but recent studies (Smith et al 2016; Williamson 2016) are taking up this crucial 
issue and raising the question of whether large-scale negative emissions are even feasible. 

100 Our aim is to contribute to understanding of the threshold-required rate of CO 2 emissions 
reduction via an approach that is transparent to non-scientists. We consider the potential for 
reductions of non-C02 GHGs to minimize the human-made climate forcing, the potential for 
improved agricultural practices to store more soil carbon, and the potential drawdown of 
atmospheric CO 2 from reforestation and afforestation. Quantitative examination reveals the 
105 merits of these actions to ameliorate demands on fossil fuel CO 2 emission phasedown, but also 
the limitations, thus clarifying the urgency of government actions to rapidly advance the 
transition to carbon-free energies to meet the climate stabilization targets they have set. 

We first describe the status of global temperature change and then summarize the principal 
climate forcings that drive long tenn climate change. We show that observed global warming is 
110 consistent with knowledge of changing climate forcings, Earth’s measured energy imbalance, 
and the canonical estimate of climate sensitivity, i.e., about 3°C 1 global warming for doubled 
atmospheric CO 2 . We illustrate updates of GHG observations and calculate a notable 
acceleration during the past decade of the growth rate of GHG climate forcing. For future fossil 
fuel emissions we consider both the IPCC Representative Concentration Pathways (RCP) 

115 scenarios, and simple emission growth rates that are helpful for determination of the plausibility 
of required emission changes. We use a precisely defined Green’s function calculation of global 
temperature with canonical climate sensitivity for each emissions scenario, thus allowing us to 
determine the amount of CO 2 that must be extracted from the air - effectively the climate debt - 
to achieve the targets of returning atmospheric CO 2 to less than 350 ppm or limiting global 
120 warming to less than 1.5°C above preindustrial levels. We discuss alternative extraction 

technologies and their estimated costs, and finally we consider the potential alleviation of CO 2 
extraction requirements that might be obtained via special efforts to reduce non-CCF GHGs. 


1 IPCC (2013) finds that 2xC0 2 equilibrium sensitivity is likely in the range 3 ± 1.5°C, as was estimated by Charney 
et al. (1979). Median sensitivity in recent model inter-comparisons is 3.2°C (Andrews et al 2012; Vial et al 2013). 
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Fig. 2. Global surface temperature relative to 1880-1920 based on GISTEMP analysis (Appendix A), (a) 
Annual and 5-year means since 1880, (b) 12- and 132-month running means since 1970. Black squares in 
(b) are calendar year (Jan-Dee) year means used to construct (a), (b) uses data through August 2016. 


2. Global Temperature Change 

The United Nations 1992 Framework Convention on Climate Change (UNFCCC 1992) stated its 
objective as ‘ . . .stabilization of GFIG concentrations in the atmosphere at a level that would 
prevent dangerous anthropogenic interference with the climate system’. The 15 th Conference of 
the Parties (Copenhagen Accord 2009) concluded that this objective required a goal to ‘ . . .reduce 
global emissions so as to hold the increase of global temperature below 2°C. . . ’ and that the 2015 
Conference of the Parties should consider the possibility of strengthening the temperature limit 
to below 1.5°C. Indeed, the Paris Agreement (2015) modified the objective ‘to holding the 
increase of global average temperature to well below 2°C above pre- industrial levels and to 
pursue efforts to limit the temperature increase to 1.5°C above preindustrial levels...’. 

Defining a target for limiting human interference with natural climate requires quantitative 
assessment of ongoing and paleo temperature changes, with the latter especially helpful for 
characterizing long-term ice sheet and sea level response versus temperature. We examine the 
modern period with near-global instrumental temperature data in the context of the current and 
previous (Holocene and Eemian) interglacial periods for which less precise proxy-based 
temperatures have recently emerged. The Holocene, now over 1 1,700 years in duration, has had 
relatively stable climate. The Eemian, lasting from about 130,000 to 1 15,000 years ago, was 
moderately warmer than the Holocene. 


2.1. Modern Temperature 

The several analyses of temperature change since 1880 are in close agreement (Hartmann et al 
2013). Thus we can use the current GISTEMP analysis (see Supporting Information), which is 
updated monthly and available (http://www.columbia.edu/~mhsl 19/Temperature/). 

The popular measure of global temperature is the annual-mean global-mean value (Fig. 2a), 
which is publicized at the end of each year. However, as discussed by Hansen et al. (2010), the 
12-month running mean global temperature is more informative and removes monthly “noise” 
from the record just as well as the calendar year average. For example, the 12-month running 
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155 mean for the past 35 years (Fig. 2b) defines clearly the super-El Ninos of 1997-98 and 2015-16 
and the 3-year cooling after the Mount Pinatubo volcanic eruption in the early 1990s. 

Global temperature in each month for the past year has been at or near a record for the 
month. Perhaps helped by a popular “spiral” temperature visualization (Hope 2016), this has 
tended to create a popular impression that global temperature may be spiraling out of control. 

160 This series of monthly records is likely to tenninate soon and the 12-month running mean is 

expected to decline as it has after prior El Ninos. However, the year-to-date temperature is so far 
above the prior record already that even the steepest post-El Nino decline cannot prevent a 2016 
annual record temperature. 

One effect of the recent warming is to remove unequivocally the illusion of a global 
165 warming hiatus after the 1997-98 El Nino. Several studies, including Trenberth and Fasullo 
(2013), England et al. (2014), Dai et al. (2015) and Rajaratnam et al (2015), showed that 
temporary plateaus are consistent with expected long-term warming due to increasing 
atmospheric GHGs. Other analyses of this specific plateau help illuminate the roles of unforced 
climate variability and natural and human-caused climate forcings in observed climate change, 
170 with the Interdecadal Pacific Oscillation (a recurring pattern of ocean-atmosphere climate 

variability) playing a major role in the wanning slowdown (Kosaka and Xie, 2013; Meehl et al, 
2014; Fyfeet al, 2016). 

Global temperature defined by the linear fit over recent decades has now reached +1.06°C 
relative to the 1880-1920 average (Fig. 2), the 12-month running mean temperature through 
175 August 2016 is 1.30°C, and the 2016 temperature likely will be near +1.25°C. The present 
global wanning rate, based on a linear fit through the past 45 years (dashed line in Fig. 2b) is 
+0.1 8°C per decade. At this rate, the trend line of global temperature, which is a relevant 
measure of mean temperature, will reach +1.5°C in about 2040 and +2°C in the late 2060s. 
However, the wanning rate can accelerate or decelerate, depending on policies that affect GHG 
180 emissions, developing climate feedbacks, and other factors discussed below. 

2.2. Temperature during current and prior interglacial periods 

Holocene temperature has been reconstructed at centennial-scale resolution from 73 globally 
distributed proxy temperature records by Marcott et al (2013). This record shows a decline of 
185 0.6°C from early Holocene maximum temperature to a “Little Ice Age” minimum in the early 

1800s [that minimum being better defined by higher resolution data of Abram et al (2016)]. 

Concatenation of the modern and Holocene temperature records (Fig. 3) assumes, based on 
Abram et al (2016), that the 1880-1920 mean temperature is 0.1 °C wanner than the Little Ice 
Age minimum. The early Holocene maximum in the Marcott et al (2013) data is thus at +0.5°C 
190 relative to the 1880-1920 mean of the modern data. However, model simulations suggest that 
the reconstructed early Holocene maximum may be ezaggerated due to limitations of the proxy 
data, especially potential seasonality bias, as discussed by Marcott et al (2013) and Liu et al 
(2014). 

Even though we cannot be certain that the current year is wanner than any single year earlier 
195 in the Holocene due to centennial smoothing of the Holocene stack and original resolution of the 
underlying proxy records (Marcott et al 2013), we conclude that the ongoing global wanning 
trend (1.06°C over 115 years, Fig 2b) is already well above prior centennially smoothed 
Holocene temperature. Further, we suggest that these smoothed temperatures are relevant to 
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Fig. 3. Estimated average global temperature for the last interglacial (Eemian) period (McKay et al 20 1 1 ; 
Clark and Huybers 2009; Turney and Jones 2010), the centennially-smoothed Holocene (Marcott et al 
2013) temperature as a function of time, and the 1 1 -year mean of modem data (Fig. 2). Vertical 
downward arrows indicate likely overestimates (see text). 

important climatic features that change on long time scales, such as ocean warming (von 
Schuckmann et al 2016), ice sheet stability (DeConto and Pollard 2016), shifting climatic zones 
(Seidel et al 2008), and the frequency of climate extremes (Hansen and Sato 2016). The formal 
2a (95% confidence) uncertainty in the Marcott et al (2013) Holocene temperature curve is only 
~0.25°C. Although total uncertainty is larger, because of issues such as discussed by Liu et al 
(2014), those uncertainties tend to push the early Holocene temperatures lower, increasing the 
gap between today’s temperature and early Holocene temperature (Marcott and Shakun 2015). 

We also conclude that the modern trend line of global temperature crossed the early 
Holocene (smoothed) temperature maximum (+0.5°C) already in about 1985. This conclusion 
receives support from the accelerating rate of sea level rise, which approached a rate of 3 
mm/year at about that date (Fig. 29 of Hansen et al 2016 shows a relevant concatenation of 
measurements). Such a high rate of sea level rise, which equates to 3 meters per millennium, far 
exceeds rates of Holocene sea level rise except in the earliest Holocene when melt was still 
coming from the final decay of mid-latitude ice sheets (Dutton et al 2015). 

The Framework Convention (UNFCCC 1992) and Paris Agreement (2015) define goals 
relevant to ‘preindustrial’ temperature, but do not define that period. We use 1880-1920, the 
earliest time with good global coverage of instrumental data, as the zero-point for temperature 
anomalies. Alternatively, one might argue for defining preindustrial as the Little Ice Age 
minimum temperature, but the deep ocean did not have time to reach equilibrium with those brief 
conditions, and global mean Little Ice Age temperature was probably only -0.1C cooler than the 
1880-1920 mean (Abram et al 2016). 

The important point is that the relevant mean global temperature has already risen out of the 
centennial Holocene range. Global warming is already having substantial adverse climate 
impacts (IPCC 2014), including extreme events (NAS 2016), and there is widespread agreement 
that 2°C warming would commit the world to multi-meter sea level rise (Levermann et al 2013; 
Clark et al 2016), and a case has been made that this could unfold within 50-150 years (Hansen 
et al 2016). 

The prior interglacial period, the Eemian, was warmer than the Holocene and sea level 
reached heights 6-9 m (20-30 feet) higher than today (Dutton et al (2015). McKay et al (2011) 
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estimated peak Eemian annual global ocean SST as +0.7°C ± 0.6°C, while models, as described 
235 by Masson-Delmotte (2013), give more confidence to the lower part of that range. Global ocean 
SST response to climate forcings is typically 70-75% as large as the global mean (land + ocean) 
surface temperature response and that same proportion is found empirically in the warming of 
the past century (http://www.columbia.edu/~mhsl 1 9/Tempera ture/T_moreFigs/). Thus the 
McKay et al data are equivalent to a global Eemian temperature +1°C relative to the Holocene. 
240 Clark and Huybers (2009) and Turney and Jones (2010) estimated global temperature in the 

Eemian as 1.5-2°C warmer than the Holocene (Fig. 3), but Bakker and Renssen (2014) analyzed 
the likely error in Eemian temperature estimates caused by the assumption that maximum 
Eemian temperatures at all proxy temperature sites occurred simultaneously and also the effect 
of proxy biases towards summer conditions, concluding that these biases could exaggerate 
245 Eemian temperature by 1 . 1 ± 0.4°C. Thus, consistent with the discussion of Masson-Delmotte 
(2013), we conclude that mean Eemian temperature was probably about 1°C warmer than the 
Holocene. Given growing indications, discussed above, that the early Holocene was little 
warmer than the pre-industrial (1880-1920) period, we conclude that Eemian global temperature 
was not much more than +1°C relative to 1880-1920 global temperature. 

250 These considerations add to the question of whether 2°C, or even 1.5°C, is an appropriate 
target to protect the well-being of young people and future generations, as modeling projections 
compared to these targets usually include only fast-feedback processes. Indeed, Hansen et al 
(2008) concluded “If humanity wishes to preserve a planet similar to that on which civilization 
developed and to which life on Earth is adapted, paleoclimate evidence and ongoing climate 
255 change suggest that C0 2 will need to be reduced from its (then) current 385 ppm to at most 350 
ppm, but likely less than that.” And further “If the present overshoot of the target C0 2 is not 
brief, there is a possibility of seeding irreversible catastrophic effects.” 

A danger of the 1.5°C and 2°C temperature targets is that they are far above the Holocene 
temperature range. If such temperature levels are allowed to long exist they will spur “slow” 

260 amplifying feedbacks (Hansen et al 2013; Rohling et al 2013; Masson-Delmotte et al 2013), 

which may have potential to run out of humanity’s control. The most threatening slow feedback 
likely is ice sheet melt and consequent sea level rise, but there are other risks in pushing the 
climate system far out of its Holocene range. Methane release from melting permafrost and 
methane hydrates is also a potentially important feedback, for example, although there are large 
265 gaps in our understanding of this feedback including its time-scale (O’Connor et al 2011). 

Thus in this paper we examine the fossil fuel emission reductions required to restore 
atmospheric C0 2 to 350 ppm or less, so as to keep global temperature close to the Holocene 
range, in addition to the canonical 1.5°C and 2°C targets. Quantitative investigation requires 
consideration of Earth’s energy imbalance, changing climate forcings, and climate sensitivity. 
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Fig. 4. Estimated effective climate forcings (update of Hansen et al 2005 through 2015). Forcings are 
based on actual changes of each gas, except CH 4 -induced changes of 0 3 and stratospheric H 2 0 are 
included in the CH 4 forcing. Oscillatory and intermittent natural forcings (solar irradiance and volcanoes) 
are excluded. CFCs include not only chlorofluorocarbons, but all Montreal Protocol Trace Gases 
(MPTGs) and Other Trace Gases (OTGs). 

3. Global Climate Forcings and Earth’s Energy Imbalance 

The dominant human-caused drivers (forcings) of climate change are changes of atmospheric 
GHGs and aerosols. GHGs absorb Earth’s infrared (heat) radiation, thus serving as a “blanket” 
that warms Earth’s surface. Aerosols, fine particles in the air that cause visible air pollution, 
both reflect and absorb solar radiation, but reflection of solar energy to space is their dominant 
effect, so they cause a cooling that partly offsets GHG warming. Estimated forcings (Fig. 4), an 
update of Fig. 28b of Hansen et al (2005), are similar to those of Myhre et al (2013) in the most 
recent IPCC report (IPCC 2013). 

Climate forcings in Fig. 4 are the planetary energy imbalance caused by preindustrial-to- 
present change of each atmospheric constituent. The CH 4 forcing includes its indirect effects, as 
increasing atmospheric CH 4 causes tropospheric ozone (O3) and stratospheric water vapor to 
increase (Myhre et al 2013). Uncertainties in the forcings, discussed by Myhre et al (2013), are 
typically 10-15% for GHGs. Uncertainty in the aerosol forcing, described by a probability 
distribution function (Boucher et al 2013), is of order 50%. Our estimate of aerosol + surface 
albedo forcing (-1.2 W/m 2 ) differs from the -1.5 W/m 2 of Hansen et al (2005), as discussed 
below, but both are within the range of the distribution function of Boucher et al (2013). 

The positive net forcing (Fig. 4) causes Earth to be out of energy balance, with more energy 
coming in than going out, which drives slow global warming. Eventually Earth will become hot 
enough to radiate to space an amount of energy matching absorbed sunlight. However, because 
of the ocean’s great thermal inertia (heat capacity), full atmosphere-ocean response to the forcing 
requires a long time: atmosphere-ocean models suggest that even after 100 years only 60-75% of 
the surface warming for a given forcing has occurred, the remaining 25-40% still being “in the 
pipeline” (Hansen et al 201 1; Collins et al 2013). Moreover, we will outline in the next section 
that global warming can activate “slow” feedbacks, such as changes of ice sheets or melting of 
methane hydrates, so the time for the system to reach a fully equilibrated state is even longer. 

GHGs have been increasing for more than a century and Earth has partially warmed in 
response. Earth’s energy imbalance is the portion of the forcing that has not yet been responded 
to. This imbalance thus defines additional global warming that will occur without further change 

Heat Storage in Upper 2000 Meters of Ocean 
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Fig. 5. Ocean heat uptake in upper 2 km of ocean during 1 1 years 2005-2015 using analysis method of 
von Schuckmann and LeTraon (2011). Heat uptake in W/m 2 (0.5 and 0.7) refer to global (ocean + land) 
area, i.e., it is the contribution of the upper ocean to the heat uptake averaged over the entire planet. 

310 of forcings. Earth’s energy imbalance can be measured by monitoring ocean subsurface 

temperatures, because almost all excess energy coming into the planet goes into the ocean (von 
Schuckmann et al 2016). Most of the ocean’s heat content change occurs in the upper 2000 m 
(Levitus et al 2012), which has been well measured since 2005 when the distribution of diving 
Argo floats achieved good global coverage (von Schuckmann and Le Traon 2011). 

315 Earth’s energy imbalance was about +0.6 W/m 2 during 2005-2010 (Hansen et al 2011) as 
inferred from heat gain in the upper 2 km of ocean using the von Schuckmann and Le Traon 
(2011) analysis and adding the smaller heat gain by the deep ocean (Purkey and Johnson 2013), 
continents, atmosphere, and net melting of sea ice and land ice. Accounting for the declining 
solar irradiance during 2005-2010 (Fig. A3), Hansen et al (2011) inferred that the energy 
320 imbalance with the solar cycle effect removed was +0.75 ± 0.25 W/m 2 . Here we update the von 
Schuckmann and Le Traon analysis with data for 2005-2015 (Fig. 5) finding 0.7 W/m 2 heat 
uptake in the upper 2000 m of the ocean. The 1 1-year period now available (2005-2015) should 
practically eliminate a solar cycle influence. The small heat gains noted above add of order 0. 1 
W/m 2 to the global heat gain (Rhein et al 2013), so our current analysis is consistent with a 
325 planetary energy imbalance of +0.75 ± 0.25 W/m 2 . The value 0.75 W/m 2 is near the middle of a 
range of estimates by several investigators (von Schuckmann et al 2016; Trenberth et al 2016). 

4. Climate Sensitivity, a Consistency Check and Slow Feedbacks 

Climate sensitivity has been a fundamental issue since at least the 19 th century when Tyndall 
330 (1861) and Arrhenius (1896) stimulated interest in the effect of a CO 2 change on climate. 

Doubled atmospheric C0 2 , a forcing of about 4 W/ 111 2 , is now a standard forcing in studies of 
climate sensitivity. The Charney et al (1979) study concluded that equilibrium sensitivity, i.e., 
global wanning after a time sufficient for the planet to restore energy balance with space, was 
3°C ± 1.5°C for 2*C0 2 or 0.75°C per W/m 2 forcing. The central value found in a wide range of 
335 modern climate models (Flato et al 2013) and in empirical paleoclimate studies remains 3°C for 
2 xCC> 2, but with an uncertainty that is still of order 1°C (Rohling et al 2012a). 

An important consistency check is obtained by comparing the estimated net climate forcing 
(2.5 W/m 2 , Fig. 4), Earth’s energy imbalance (-0.75 W/m 2 ), observed global warming, and 
climate sensitivity. Observed wanning since 1880-1920 is 1.06°C with the effect of El Nino/La 
340 Nina oscillations removed (Fig. 2b). Global warming between 1700-1800 and 1880-1920 was 
~0.1°C (Abram et al 2016; Marcott et al 2013), so 1750-2015 warming was ~1.16°C. Taking 
climate sensitivity as 0.75°C per W/m 2 forcing, global warming of 1.16°C implies that 1.55 
W/m 2 of the total 2.5 W/m 2 forcing has been “used up” to cause observed warming. Thus 0.95 
W/m 2 forcing should remain to be responded to, i.e., the expected planetary energy imbalance is 
345 0.95 W/m 2 , reasonably consistent with the observed 0.75 ± 0.25 W/m 2 . If we instead use the 

aerosol + surface albedo forcing -1.5 W/m 2 estimated by Hansen et al (2005), the net climate 
forcing is 2.2 W/m 2 and the forcing not responded to is 0.65 W/m 2 , which is also within the 
observational error of Earth’s energy imbalance. 

An important matter to bear in mind is that the sensitivity 3°C for 2* CO 2 (0.75°C per W/m 2 ) 
350 is the “fast-feedback” climate sensitivity, i.e., it does not include “slow” climate feedbacks that 
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will occur if global temperature long remains above the Holocene level (Hansen et al 2008; 
Rohling et al 2012a). Slow feedbacks include large-scale shrinking of ice sheets as Earth wanns 
and the enhanced release of GHGs as the ocean, soil, and continental shelves warm. These slow 
feedbacks are strongly amplifying, indeed, they are the reason that natural long-term climate 
355 oscillations are so large in response to even small long-term global-average forcings (Rohling et 
al 2012b; Masson-Delmotte et al 2013). 

The fast-feedback climate sensitivity is the appropriate sensitivity to use in interpretation of 
recent climate change, because we use observed change of GHGs and because ice sheet change 
so far is small. However, the need to avoid the emergence of slow feedbacks motivates the 
360 criterion that energy balance should be restored at a global temperature close to and eventually 
within the Holocene range (Hansen et al. 2008, 2013). 

Earth’s present energy imbalance is causing heat to accumulate in the ocean, where it 
contributes to melting of ice shelves (Rignot et al 2013). Rising temperatures also increase the 
risk of CO 2 and CH 4 release from drying soils, thawing permafrost (Schadel et al 2016; Schuur et 
365 al 2015) and wanning continental shelves (Kvenvolden 1993, Judd et al 2002). Time scales for 
the slow feedbacks are not well established, but recent modeling and empirical evidence suggest 
that substantial ice sheet and sea level changes could occur within periods as short as several 
decades (Rohling et al 2013; Pollard et al 2015; Hansen et al 2016). If large planetary energy 
imbalance continues, there is a danger that the warming driving slow feedbacks will be so far 
370 advanced that consequences such as large sea level rise proceed out of humanity’s control. 

Quantification of requirements for stabilizing climate depends on knowledge of ongoing 
changes of the two largest GHG forcings, C0 2 and CH 4 . It is also necessary to understand how 
we are changing GHG emissions directly through industrial and agricultural activities 
(designated ‘anthropogenic emissions’ and included in the SRES and RCP scenarios) and 
375 indirectly through climate change (the slow feedbacks noted above, designated somewhat 

paradoxically ‘natural emissions’ changes, and not included in the SRES and RCP scenarios). 
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Fig. 6. (a) Global C0 2 annual growth based on NOAA data ('http://www.esrl.noaa.gov/gmd/ccgg/trends/) . 
Dashed curve is for a single station (Mauna Loa). Red curve is monthly global mean relative to the same 
month of prior year; black curve is 12-month running mean of red curve, (b) C0 2 growth rate is highly 
correlated with global temperature, the C0 2 change lagging global temperature change by 8 months. 


5. Observed CO 2 and CH 4 Growth Rates 

Annual increase of atmospheric C0 2 , averaged over a few years, grew from less than 1 ppm/year 
50 years ago to more than 2 ppm/year today (Fig. 6), with the global mean and Mauna Loa C0 2 
amounts now exceeding 400 ppm (Betts et al 2016). The large oscillations of the annual growth 
are correlated with global temperature and with the El Nino/La Nina cycle 2 . Correlations are 
calculated for the 12-month running means, which effectively remove the seasonal cycle and 
monthly noise. Maxima of the CO 2 growth rate lag global temperature maxima by ~8 months 
(Fig. 6b) and lag Nino3.4 [latitudes 5N-5S, longitudes 120-170W] temperature by ~10 months. 
These lags imply that the current C0 2 growth spike (Fig. 6 uses data through July 2016), 
associated with the 2015-16 El Nino, may not have reached its maximum yet, as Nino3.4 peaked 
in December 2015 and global temperature peaked in February 2016. 
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Fig. 7. Fossil fuel C0 2 emissions (left scale) and airborne fraction, i.e., the ratio of observed atmospheric 
C0 2 increase to fossil fuel C0 2 emissions. 


2 One mechanism for greater than normal atmospheric C0 2 growth during El Ninos is the impoverishment of 
nutrients in equatorial Pacific surface water and thus reduced biological productivity that result from reduced 
upwelling of deep water (Chavez et al., 1999). However, the El Nino/La Nina cycle seems to have an even greater 
impact on atmospheric C0 2 via the terrestrial carbon cycle through effects on the water cycle, temperature, and fire, 
as discussed in a large body of literature (referenced, e.g., by Schwalm et al., 2011). 
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Fig. 8. Global CH 4 from Dlugokencky (2016), NOAA/ESRL (www.esrl.noaa.gov/gmd/ccgg/trends_ch4/). 
End months for the three indicated slopes are January 1984, May 1992, August 2006, and January 2016 

Growth of airborne CCE (defined as the increase in atmospheric CCE above preindustrial 
levels) appears to be about half of fossil fuel emissions (Fig. 7), the remaining portion being net 
uptake by the ocean and biosphere (Ciais et al 2013). Here we use the Keeling et al. (1973) 
definition of airborne fraction, which is the ratio of quantities that are known with good 
accuracy: the annual increase of CCE in the atmosphere and the annual amount of CO 2 injected 
into the atmosphere by fossil fuel burning. The data reveal that, even as fossil fuel emissions 
have increased by a factor of four over the past half century, the ocean and biosphere have 
continued to take up about half of the emissions (Fig. 7, right-hand scale). This seemingly 
simple relation between emissions and atmospheric CO 2 growth is not predictive as it depends on 
the growth rate of emissions being maintained, and, indeed, it is not expected to continue in 
cases with major changes in the emission scenario, so we use a carbon cycle model in Section 7 
to compute atmospheric CO 2 as a function of emission scenario. 

Atmospheric CH 4 stopped growing between 1998 and 2006, indicating that its sources and 
sinks were nearly in balance, but growth resumed in the past decade (Fig. 8). Growth of CH 4 
exceeds 10 ppb/year in 2014 and 2015, almost as fast as in the 1980s). Turner et al. (2016) 
suggest that increased fossil fuel emissions in the U.S. may be a major cause of renewed global 
CH 4 growth. However, CH 4 isotope data imply that resumed growth was mainly from wetlands, 
especially in the tropics but with a contribution from high latitudes of the Northern Hemisphere 
(Bousquet et al 201 1; Dlugokencky et al 2011). The CH 4 changes over the past two decades are 
driven primarily by changes in emissions as observations of CH3CCI3 show very little change in 
the atmospheric sink for CH 4 ( Montzka et al. 2011: Holmes et al. 2013) . Future changes in the 
sink, however, are expected to lead to increased atmospheric CH 4 separate from emission 
changes, but these are difficult to project in the RCP scenarios ( Voulgarakis et al. 2013 ). 

The continued growth of atmospheric CO 2 and the reaccelerating growth of CH 4 raise 
important questions related to prospects of stabilizing climate. How consistent are scenarios for 
phasing down climate forcing with reality revealed by observational data? What changes to 
emissions are required to stabilize climate? We address these issues below. 
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Fig. 9. GHG climate forcing growth rate with historical data being 5-year running means, except data for 
2014 and 2015 are 3- and 1-year means, (a) includes scenarios used in IPCC AR3 and AR4 reports, and 
(b) has AR5 scenarios. N 2 0, MPTGs and OTGs (Montreal Protocol Trace Gases and Other Trace Gases) 
435 data are from NOAA/ESRL Global Monitoring Division. 

6. GHG Climate Forcing Growth Rates and Emission Scenarios 

Insight is obtained by comparing the growth rate of GHG climate forcing based on observed 
GHG amounts with past and present GHG scenarios. We examine forcings of IPCC SRES 
440 (2000) scenarios used in the AR3 and AR4 reports (Fig. 9a) and RCP scenarios (IPCC 2013) 

used in the AR5 report (Fig. 9b). We include the “alternative scenario” of Hansen et al (2000) in 
which CO2 and CH 4 emissions decline such that global temperature stabilizes near the end of the 
century. 3 We use the same radiation equations for observed GHG amounts and scenarios, so 
errors in the radiation calculations do not alter the comparison. Equations for GHG forcings are 
445 from Table 1 of Hansen and Sato (2004) with the CH 4 forcing using an efficacy factor 1.4 to 
include effects of CH 4 on tropospheric O3 and stratospheric H 2 0 (Hansen et al 2005). 

The growth of GHG climate forcing peaked at -0.05 W/nr/year (5 W/m 2 /century) in 1978- 
1988, then falling to a level 10-25% below IPCC SRES (2000) scenarios during the first decade 
of the 21 st century (Fig. 9a). The decline was due to (1) decline of the airborne fraction of C0 2 
450 emissions (Fig. 7), (2) slowdown of CH 4 growth (Fig. 8), and (3) the Montreal Protocol, which 
initiated phase-out of gases that destroy stratospheric ozone. 

The situation in 2000 seemed ripe for a pathway to climate stabilization more rapid than any 
of the IPCC scenarios. The slowing growth of forcings was partly good fortune, but also due to 


3This scenario is discussed by Hansen and Sato (2004). CH 4 emissions decline moderately, producing a small 
negative forcing. C0 2 emissions (not captured and sequestered) are assumed to decline until in 2100 fossil fuel 
10 emissions just balance uptake of C0 2 by the ocean and biosphere. C0 2 emissions continue to decline after 2100. 
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Fig. 10. Fossil fuel emission scenarios. Scenarios in (a) have constant emissions in 2015-2020 and then 
simple specified rates of emission increase or decrease. IPCC (2013) RCP scenarios are shown in (b). 


the prescience of the Montreal Protocol, whose design and implementation to reduce the ozone- 
depleting gases also allowed it to be used to slow or reverse growth of some other GFIGs. The 
alternative scenario aimed to extend the downward trend in the growth rate of climate forcing by: 
(1) slowing the growth of CO 2 emissions, as may occur with a substantial rising price on carbon 
emissions to accelerate development of carbon-free energies, (2) a global effort to reduce CH 4 
emissions, (3) continued use and tightening of the Montreal Protocol to constrain trace GFIGs. 
The slowly decreasing forcing of this alternative scenario would have kept global wanning well 
below 1.5°C, for climate sensitivity 0.75°C/W/m 2 (Fig. A4). 

However, in reality, in the absence of a universally rising carbon price and substantial 
support for energy research and development, global fossil fuel C0 2 emissions accelerated, from 
1.5%/year in 1973-2000 to ~ 2 . 5%/year after 2000 (Figs. 1 and SI). The growth rate of GHG 
forcing now exceeds the alternative scenario by ~70% (Fig. 9a). New scenarios must begin from 
current reality, and, as a consequence or recent growth, ambitious targets for limiting global 
warming now require much steeper emissions reductions. 

The new IPCC (2013) RCP scenarios (Fig. 9b) initiate in 201 1 and fan out into an array of 
potential futures driven by assumptions about energy demand, fossil fuel prices, and climate 
policy, chosen to be representative of an extensive literature on possible emissions trajectories 
(Moss et al 2010; van Vuuren et al 201 1; Meinshausen et al 2011). Numbers on the RCP 
scenarios (8.5, 6.0, 4.5 and 2.6) refer to the GHG climate forcing (W/m 2 ) in 2100. 

As a complement to RCP scenarios, we define scenarios simply by percent annual emission 
decrease or increase. We consider rates -6%/year, -3%/year, constant emissions, and +2%/year; 
emissions stop increasing in the +2%/year case when they reach 25 Gt/year (Fig. 10a). Scenarios 
with decreasing emissions are preceded by constant emissions for 2015-2020, in recognition that 
some time is required to achieve policy change and implementation. Note similarity of RCP 2.6 
with -3%/year, RCP 4.5 with constant emissions, and RCP 8.5 with +2%/year (Fig. 10). 

Scenario RCP2.6 has the world moving into negative growth of GHG forcing 25 years from 
now (Fig. 9b), through rapid reduction of GHG emissions and C0 2 capture and storage. Already 
in 2015 there is a huge gap between reality and RCP2.6. Closing the gap (0.01 W/m 2 ) between 
actual growth of GHG climate forcing in 2015 and RCP2.6 (Fig. 9b), with C0 2 alone, would 
require extraction from the air of more than 0.7 ppm of CO 2 or 1.5 GtC in the single year (2015). 
We discuss the plausibility and estimated costs of scenarios with CO 2 extraction in Section 9. 
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Fig. 11 . (a) Atmospheric C0 2 for emission scenarios of Fig. 10a. (b) Atmospheric C0 2 including effect 
of C0 2 extraction that increases linearly after 2020 (after 2015 in +2%/year case). 1 ppm is -2.12 GtC. 


7. Future CO 2 for Assumed Emission Scenarios 

We must model Earth’s carbon cycle, including ocean uptake of carbon, deforestation, forest 
regrowth and carbon storage in the soil, for the purpose of simulating future atmospheric C0 2 as 
a function of fossil fuel emission scenario. Fortunately, the convenient dynamic-sink pulse- 
response function version of the well-tested Bern carbon cycle model (Joos et al 1996) does a 
good job of approximating more detailed models, and it produces a good match to observed 
industrial-era atmospheric C0 2 . Thus we use this relatively simple model, described elsewhere 
(Joos et al 1996; Kharecha and Hansen 2008 and references therein), to examine the effect of 
alternative fossil fuel use scenarios on the growth or decline of atmospheric C0 2 . For land use 
C0 2 emissions in the historical period, we use the values labeled Houghton/2 by Hansen et al 
(2008), which were shown in the latter publication to yield good agreement with observed C0 2 . 
We use fossil fuel C0 2 emissions data for 1850-2013 from Boden et al (2016). BP fuel 
consumption data for 2013-2015 is used with the fractional annual changes for each nation to 
allow extension of the Boden analysis through 2015. Emissions were almost flat from 2014 to 
2015, due to economic slowdown and increased use of low-carbon energies, but, even if a peak 
in global emissions is near, substantial decline of emissions is dependent on acceleration in the 
transfonnation of energy production and use (Jackson et al 2016).. 

The scenarios shown in Figs. 10a and 11a are the baseline cases without any anthropogenic 
C0 2 removal. We illustrate five cases with C0 2 removal in Fig. lib that achieve atmospheric 
C0 2 targets of either 350 ppm or 450 ppm in 2100, with cumulative removal amounts listed in 
parentheses. The rate of C0 2 extraction in all cases increases linearly from zero in 2010 to the 
value in 2100 that achieves the atmospheric C0 2 target (350 ppm or 450 ppm). The amount of 
C0 2 that must be extracted from the system exceeds the difference between the atmospheric 
amount without extraction and the target amount, e.g., constant C0 2 emissions and no extraction 
yields 546 ppm for atmospheric C0 2 in 2100, but to achieve a target of 350 ppm the required 
extraction is 328 ppm, not 546 - 350 = 196 ppm. The well-known reason (Cao and Caldeira 
2010) is that ocean out gassing increases, and vegetation productivity and ocean C0 2 uptake 
decrease with decreasing atmospheric C0 2 , as explored in a wide range of Earth System models 
(Jones et al 2016). 
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8. Simulations of Global Temperature Change 

Analysis of future climate change and policy options to alter that change must address various 
525 uncertainties. One useful way to treat uncertainty is to use results of many models and construct 
probability distributions (Collins et al 2013). Such distributions have been used to estimate the 
remaining budget for fossil fuel emissions for a specified likelihood of staying under a given 
global wanning limit and to compare alternative policies for limiting climate forcing and global 
warming (Rogelj et al 2016a,b). 

530 Our aim here is a fundamental, transparent calculation that clarifies how future warming 

depends on the rate of fossil fuel emissions. We use best estimates for fundamental uncertain 
quantities such as climate sensitivity. If these estimates are accurate, actual temperature should 
have about equal chances of falling higher or lower than the calculated value. Among the 
important uncertainties in projections of future climate forcings and climate change are climate 
535 sensitivity, the effects of ocean mixing and dynamics on the climate response function discussed 
below, and aerosol climate forcing. We provide all defining data so that others can easily repeat 
calculations with alternative choices. 

We calculate global temperature change T at time t in response to any climate forcing 
scenario using the Green’s function (Hansen 2008) 

540 

T(t) = JR(t) [dF/dt] dt (1) 

where R(t) is the product of equilibrium global climate sensitivity and the dimensionless climate 
response function (percent of equilibrium response), dF/dt is the annual increment of net forcing, 
545 and the integration begins before human-made climate forcing is substantial. Our response 
function reaches 75% response in 100 years, a rate that Hansen et al (2011) conclude is 
representative of the real world, based on observations of Earth’s energy imbalance; this 
imbalance is an immediate consequence of the time required for the ocean surface temperature to 
respond to changing climate forcing. Our results can be exactly reproduced, or altered with 
550 alternative choices for climate forcings, climate sensitivity and response function, as we tabulate 
the forcings in Table SI and the response function is exactly defined. 4 

We use equilibrium fast-feedback climate sensitivity % °C per W/m 2 (3°C for 2x002). This 
is consistent with current climate models (Collins et al 2013: Flato et al 2013) and paleoclimate 
evidence (Rohling et al 2012a; Masson-Dehnotte et al 2013; Bindoff and Stott 2013). 

555 CO 2 is the dominant forcing in scenarios for future climate. The growth of non-C02 GHG 
climate forcing is likely to be even smaller, relative to CO 2 forcing, than it has been in recent 
decades (Fig. 9), especially if there is a strong effort to limit climate change. Indeed, recent 
agreement to use the Montreal Protocol (2016) to phase down emissions of minor trace gases 
should cause added forcing of Montreal Protocol Trace Gases (MPTGs) + Other Trace Gases 
560 (OTGs) (red region in Fig. 9) to become near zero or slightly negative, thus at least partially off- 
setting growth of the N 2 0 climate forcing. Some N 2 0 increase may be inevitable, because its 
emissions are largely associated with food production, and population is not expected to stabilize 
before mid-century at the earliest (Ciais et al 2013; Kroeze and Bouwman 2011). 


4We use the "intermediate" response function in Fig. 5 of Hansen et al. (2011), which gives best agreement with 
Earth's energy imbalance. Fractional response is 0.15, 0.55, 0.75 and 1 at years 1, 10, 100 and 2000 with these 
values connected linearly in log (year), cf. Fig. 5 of Hansen et al (2011). 
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Fig. 12. Climate forcings used in our climate simulations; Fe is effective forcing, as discussed in 
connection with Fig. 4. (a) Future GF1G forcing uses four alternative fossil fuel emission growth rates. 

(b) GF1G forcings are altered based on CO 2 extractions of Fig. 11. 

The net effect of nitrogen emissions is complex because of both diminishing and amplifying 
feedbacks (Kroeze and Bouwman 2011), e.g., fertilizers can increase uptake of carbon by the 
biosphere and affect tropospheric O3, but it is expected that more efficient use of fertilizers can 
reduce emissions and N 2 0 growth (Liu and Zhang 2011). CH 4 is responsible for the largest non- 
CO 2 GHG forcing, with potential to significantly exacerbate or alleviate the magnitude of global 
warming, so we address the range of CH 4 possibilities in Section 11. Here we use RCP6.0 for 
the non-C0 2 GHGs, a scenario in which warming by these gases, compared to C0 2 , is small. 

We take tropospheric aerosol plus surface albedo forcing as -1.2 W/m 2 in 2015, presuming 
the aerosol and albedo contributions to be -1 W/m 2 and -0.2 W/m 2 , respectively. We assume a 
small increase this century as global population rises and increasing aerosol emission controls in 
emerging economies tend to be offset by increasing development elsewhere, so aerosol + surface 
forcing is -1.5 W/nr in 2100. The temporal shape of the historic aerosol forcing curve (Table 
SI) is from Hansen et al (2011), which in turn was based on the Novakov et al (2003) analysis of 
how aerosol emissions have changed with technology change. 

Historic stratospheric aerosol data (Table SI, annual version), an update of Sato et al (1993), 
include moderate 21 st century aerosol amounts (Bourassa et al 2012). Future aerosols, for 
realistic variability, include three volcanic eruptions in the rest of this century with properties of 
the historic Agung, El Chichon and Pinatubo eruptions, and a background stratospheric aerosol 
forcing -0.1 W/m 2 . This leads to mean stratospheric aerosol climate forcing -0.25 W/m 2 for the 
21 st century, similar to the prior century. Reconstruction of historical solar forcing (Coddington 
et al 2015; Kopp et al 2016), based on data in Fig. A3, is extended with an 1 1-year cycle. 
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Fig. 13. Simulated global temperature for forcings of Fig. 12. Observations as in Fig. 2. Gray area is 2a 
(95% confidence) range for centennially-smoothed Flolocene maximum, but there is further uncertainty 
about the magnitude of the Flolocene maximum, as noted in the text and discussed by Liu et al (2014). 

Individual and net climate forcings for the several fossil fuel emission reduction rates are 
shown in Fig. 12a, c. Scenarios with linearly growing CCL extraction at rates required to yield 
350 or 450 ppm airborne CO 2 in 2100 are in Fig. 12b, d. These forcings and the assumed climate 
response function define expected global temperature for the entire industrial era (Fig. 13). 

A stark summary of alternative futures emerges from Fig. 13. If emissions grow 2%/ year, 
modestly slower than the 2.6%/year growth of 2000-2015, warming reaches ~3°C by 2100. 
Warming is close to 2°C if emissions are constant until 2100. Furthermore, both scenarios 
launch Earth onto a course of more dramatic change well beyond the initial 2-3°C global 
warming, because: (1) warming continues beyond 2100 as the planet is still far from equilibrium 
with the climate forcing, and (2) warming of 2-3°C would unleash strong slow feedbacks, 
including melting of ice sheets and increases of GHGs, thus continuing growing climate change. 
Reducing global emissions at a rate of 3 %/year (or more steeply) maintains global warming at 
less than 1.5°C above preindustrial, but the temperature at the end of the century continues to be 
0.5°C or more above the prior Holocene maximum with consequences that are difficult to 
foresee, especially due to the likelihood of initiating substantial amplifying slow feedbacks. 

Desire to avoid slow feedbacks, including ice sheet shrinkage and sea level rise, spurs the 
need to get global temperature back into the Holocene range. This goal needs to be achieved on 
the time scale of a century or less, as paleoclimate evidence indicates that the response time of 
sea level to climate change is 1-4 centuries (Grant et al 2012, 2014) for natural climate change, 
and it is unlikely that the response would be slower to a stronger, more rapid human-made 
climate forcing. The scenarios that reduce CO 2 to 350 ppm succeed in getting temperature back 
close to the Holocene maximum by 2100 (Fig. 13b), but they require extractions of atmospheric 
CO 2 that range from 72 ppm in the scenario with 6%/year emission reductions to 768 ppm in the 
scenario with +2%/year emission growth. 

Scenarios ranging from constant emissions to +2%/year emissions growth can be made to 
yield 450 ppm in 2100 via extraction of 160-600 ppm of CO 2 from the atmosphere (Fig. 12b). 
However, these scenarios still yield warming more than 1.5°C above the preindustrial level 
(more than 1°C above the early Holocene maximum). Consequences of such warming and the 
plausibility of extracting such huge amounts of atmospheric C0 2 are considered below. 
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9. CO 2 Extraction: Plausibility and Cost 

625 The above calculations show the need for extraction of CO 2 from the air, also called negative 
emissions, in addition to reducing emissions of GHGs. A goal of 100 GtC (47 ppm C0 2 ) 
extraction in the 21 st century was chosen by Hansen et al (2013), because it is comparable to net 
emissions from historic deforestation and land use (Ciais et al 2013), and thus it is likely to be 
about as much as can be achieved via relatively natural reforestation and afforestation (Canadell 
630 and Raupach 2008) and improved agricultural practices that increase soil carbon (Smith 2016). 

We differentiate between the limited carbon that can be extracted from the air by improved 
agricultural and forestry practices and additional “technological extraction” by intensive negative 
emission technologies that might be used to remediate overshoot of the C0 2 level needed to 
assure an acceptable long-term climate state. We assume that improved practices will aim at 
635 optimizing agricultural and forest carbon uptake via relatively natural approaches, compatible 
with delivering a range of ecosystem services from the land (Smith 2016; Smith et al 2016) In 
contrast, proposed technological extraction and storage of CO 2 does not have co-benefits and 
remains unproven at relevant scales (NRC 2015). Improved practices have local benefits in 
agricultural yields and forest products and services (Smith et al 2016), which may help minimize 
640 net costs. Developed countries recognize a financial obligation to less developed countries that 
have done little to cause climate change (Paris Agreement 20 15) 5 . We suggest that at least part 
of developed country support should be channeled through an agricultural and forestry program, 
with continual evaluation and adjustment to reward and encourage progress (Bustamante et al 
2014). Non-C0 2 GHGs could be included in the improved practices program. We do not 
645 estimate the program cost, but we assume that such a program will be carried out, if there is to be 
hope of stabilizing climate. Thus the costs we estimate for additional technological extraction of 
CO 2 are a minimum cost. 

Here we first reexamine the question of whether a concerted global effort on carbon storage 
in forests and soil might have potential to provide a carbon sink substantially larger than 100 GtC 
650 this century. Smith et al. (2016) estimate that reforestation and afforestation together have 
carbon storage potential of about 1 . 1 GtC/year. However, as forests mature, their uptake of 
atmospheric carbon decreases (termed “sink saturation”), thereby limiting CO 2 drawdown. 

Taking 50 years as the average time for tropical, temperate and boreal trees to experience sink 
saturation yields 55 GtC as the potential storage in forests this century. 

655 Smith (2016) shows that soil carbon sequestration and soil amendment with biochar 

compare favorably with other negative emission technologies with less impact on land use, water 
use, nutrients, surface albedo, and energy requirements, but understanding of and literature on 
biochar are limited (NRC 2015). Smith concludes that soil carbon sequestration has potential to 
store 0.7 GtC/year. However, as with carbon storage in forest, there is a saturation effect. A 
660 commonly used 20-year saturation time (IPCC 2006) would yield storage of 14 GtC soil carbon 
storage, while an optimistic 50-year saturation time would yield 35 GtC. Use of biochar to 
improve soil fertility provides additional carbon storage with potential rate as high as 0.7- 1.8 
GtC/year (Woolf et al 2010; Smith 2016). Larger industrial-scale biochar carbon storage is 
conceivable, but belongs in the category of intensive negative emission technologies, discussed 

5 Another conceivable source of financial support for C0 2 drawdown might be legal settlements with fossil fuel 
15 companies, analogous to penalties that courts have imposed on tobacco companies, but with the funds directed to 
the international "improved practices" program. 
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665 below, whose environmental impacts and costs require scrutiny. We conclude that 100 GtC is an 
appropriate estimate for potential carbon extraction via an ambitious concerted global-scale 
effort to improve agricultural and forestry practices with carbon drawdown as a prime objective. 

Copious CO 2 extraction is conceivable via other intensive negative emission technologies, 
including (1) burning of biofuels in power plants with capture and sequestration of resulting C0 2 
670 (Creutzig et al 2015), and (2) direct air capture of C0 2 and sequestration (Keith 2009; NRC 
2015), and (3) grinding and spreading of minerals such as olivine to enhance the geological 
weathering process (Taylor et al 2016). However, energy, land and water requirements of these 
technologies impose economic and biophysical limits on CO 2 extraction (Smith et al 2016). 

The popular concept of bioenergy with carbon capture and storage (BECCS) requires large 
675 areas, high fertilizer and water use, and may compete with other vital land use such as agriculture 
(Smith, 2016). Costs estimates are ~$ 1 50-350/tC for crop-based BECCS (Smith et al 2016). 

Direct air capture has less area and water needs than BECCS and no fertilizer requirement, 
but it has high energy use, has not been demonstrated at scale, and cost estimates exceed those of 
BECCS (Socolow et al 201 1; Smith et al 2016). Keith et al (2006) have argued that, with strong 
680 research and development support and industrial-scale pilot projects sustained over decades, it 
may be possible to achieve costs ~$200/tC, thus comparable to BECCS costs; however other 
assessments are higher, reaching $1400-3700/tC (NRC 2015). Carbon capture and storage 
(CCS) from a stream of nearly 100 percent CO 2 at fossil fuel burning sites is more efficient and 
thus less expensive than direct air capture, but CCS at power plants is properly included in our 
685 scenarios as one of the mechanisms competing to achieve phase-down of fossil fuel emissions, 
along with energy efficiency, renewable energies, and nuclear power. 

Enhanced weathering via soil amendment with crushed silicate rock is a candidate negative 
emission technology that also limits coastal ocean acidification as chemical products liberated by 
weathering increase land-ocean alkalinity flux (Kohler et al 2010; Taylor et al 2016). If two- 
690 thirds of global croplands were amended with basalt dust, as much as 2-5 GtC/year might be 
extracted, depending on application rate (Taylor et al 2016), but energy costs from mining, 
grinding and spreading likely reduce this by 10-25% (Moosdorf et al 2014). Although such 
large-scale enhanced weathering is speculative, there are potential co-benefits for temperate and 
tropical agroecosystems that could affect its practicality, and may put some enhanced weathering 
695 into the category of improved agricultural and forestry practices. Benefits include fertilizing of 
crops that increases yield and reduces use and cost of other fertilizers, increasing crop protection 
from insect herbivores and pathogens thus decreasing pesticide use and cost, neutralizing soil 
acidification to improve yield, and suppression of GHG (N 2 0 and C0 2 ) emissions from soils 
(Edwards et al 2016; Kantola et al 2016). Cost of enhanced weathering might be reduced by 
700 deployment with reforestation and afforestation and with crops used for BECCS. 

For cost estimates, we first consider restoration of airborne CO 2 to 350 ppm in 2100 (Fig. 

1 lb), which would keep global warming below 1.5°C and bring global temperature back close to 
the Holocene maximum by end-of-century (Fig. 13b). This scenario keeps the temperature 
excursion above the Holocene level small enough and brief enough that it has the best chance of 
705 avoiding ice sheet instabilities and multi-meter sea level rise (Hansen et al 2016). If fossil fuel 
emission phasedown of 6%/year had begun in 2013, as proposed by Hansen et al (2013), this 
scenario would have been achieved via the 100 GtC carbon extraction from improved 
agricultural and forestry practices. 
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Now, with assumption that global emissions will be comparable to today’s level through 
710 2020, Figs. 1 lb and 13b show that 6%/year emissions reduction starting in 2021 leaves a 

requirement to extract 72 ppm CO 2 (153 GtC) from the air during this century. Emission 
reductions of 3%/year leave a requirement of extracting 1 12 ppm C0 2 (Fig. 13b) by 2100. 
Constant emissions and +2%/year emissions growth would require extractions of 328 and 768 
ppm CO 2 to reach 350 ppm in 2100. 

715 The lowest cost is for the case of 6%/year emissions reduction. We assume that 100 GtC 

will be stored in the biosphere via improved agricultural and forestry practices. We do not mean 
to diminish the magnitude or cost of this task, but we must assume that it will occur if climate 
change impacts are to be minimized, and further we expect that developed countries will 
recognize their obligations to provide assistance required to achieve success. 

720 The remaining 53 GtC, at the rate $ 1 50-350/tC estimated for BECCS and other intensive 
negative emission technologies (Fig. 3f of Smith et al 2016), would cost $8-18.5 trillion, thus 
$100-230 billion per year if spread uniformly over 80 years. In contrast, continued high 
emissions, say between constant emissions and +2%/year, require extraction of 695-1628 GtC, 
which corresponds to $104-570 trillion dollars or $1.3-7 trillion dollars per year over 80 years. 6 
725 Such extraordinary cost, along with the land area, fertilizer and water requirements (Smith et al 
2016) suggest that, rather than the world being able to buy its way out of climate change, 
continued high emissions may force humanity to largely live with the climatic consequences. 

10. Climate Forcing Contribution of N011-CO2 GHGs 

730 GHG climate forcing is surging, not declining, the annual rate having increased more than 20% 
in just the past five years (Fig. A5). This recent surge in the growth rate of the GHG climate 
forcing is led by increasing growth of CH 4 , but C0 2 is by far the largest cause of continued 
growth of the GHG climate forcing (Fig. 9). Given the difficulty and cost of reducing C0 2 , we 
must ask about the alternative of reducing non-CCk GHGs. Could realistic reductions of these 
735 other gases substantially alter the CO 2 abundance required to meet a target climate forcing? 

Methane (CH 4 ) is the largest climate forcing other than CO 2 (Fig. 3). The CH 4 
atmospheric /lifetime is only about 10 years (Prather et al 2012), so there is potential to reduce 
this climate forcing rapidly if CH 4 sources are reduced. Our climate simulations, employing 
RCP6.0 for non-C02 gases, make an optimistic assumption that future CH 4 , after a moderate 
740 increase in the next few decades, will decrease from its present -1800 ppb to 1650 ppb in 2100, 
yielding a forcing -0.1 W/ 111 2 . RCP2.6 makes a more optimistic assumption: that CH 4 will 
decline monotonically to 1250 ppb in 2100, yielding a forcing -0.3 W/m 2 (relative to today’s 
1800 ppb CH 4 ), based on radiation equations identified in section 6. 


6 For reference, the United Nations global peacekeeping budget is about $10B/year. National military budgets are 
larger: the 2015 USA military budget was $596B and the global military budget was $1.77 trillion (SIPRI 2016). 
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(a) CH 4 Abundance (ppb) (b) N 2 0 Abundance (ppb) 



Fig. 14. Comparison of observed CH 4 and N 2 0 amounts and RCP scenarios. RCP 6.0 and 4.5 scenarios 
for N 2 0 overlap. Observations are from NOAA/ESRL Global Monitoring Division. 

Actual atmospheric CH 4 abundance (Fig. 14) is diverging on the high side from these 
optimistic scenarios. The downward offset (~20 ppb) of CH 4 scenarios relative to observations 
(Fig. 14) is due to the fact that RCP scenarios did not include a data adjustment that was made in 
2005 to match a revised CH 4 standard scale (E. Dlugokencky, priv comm). In addition, observed 
CH 4 is increasing more rapidly than in most scenarios. 

Carbon isotopes provide a valuable constraint on which CH 4 sources 7 contribute to the CH 4 
growth resurgence in the past decade (Fig. 7). Specifically, Schaefer et al (2016) conclude that 
the growth was primarily biogenic, thus not fossil fuel, and located outside the tropics, most 
likely ruminants and rice agriculture. Such an increasing biogenic source is consistent with 
effects of increasing population and dietary changes (Tilman and Clark 2014). These sources 
potentially could be mitigated [by changing rice growing methods (Epule et al 201 1) and 
inoculating ruminants (Eckard et al, 2010; Beil 2015)], but that would require widespread 
adoption of new technologies at the farmer level. Concerning fossil fuels, it is feasible to reduce 
CH 4 leaks, yet enhanced shale gas extraction of CH 4 may yield even greater leakage (Caulton et 
al., 2014; Petron et al., 2014; Howarth, 2015). 

Slow climate feedbacks could increase CH 4 because natural emissions from methane 
hydrates, permafrost and natural wetlands 8 are expected to increase in response to global 
warming (O’Connor et al 2010). On the other hand, as yet there is little evidence for substantial 
emissions from hydrates or permafrost (Warwick et al., 2016). Predicting such emissions has 
large uncertainty, because drought conditions eliminate wetland CH 4 emissions (while greatly 
increasing C0 2 release from soil carbon), and, in addition, CH 4 created in anoxic zones is mostly 
oxidized in the water column before reaching the atmosphere (Reeburgh, 2007). 

All reasonable effort to reduce methane is appropriate, recognizing that its mitigation effort 
will be different than that for fossil fuel C0 2 in that it should include a focus on agriculture in 
developing countries with adoption of new practices and technology at the farm level. However, 

7 Estimated human-caused CH 4 sources (Ciais et al., 2013) are: fossil fuels (29%), biomass/biofuels (11%), Waste 
and landfill (23%), ruminants (27%) and rice (11%) 

8 Wetlands compose a majority of natural CH 4 emissions and are estimated to be equivalent to about 36% of the 
anthropogenic source (Ciais et at, 2013) 
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given increasing global population and global wanning “in the pipeline,” there is an underlying 
tendency for greater emissions. The current CH 4 increases (Fig. 14a) show that the mitigation 
775 paths envisaged with RCP scenarios projecting methane decreases are not close to present 

reality. Nevertheless, it is plausible that human-caused emissions could achieve a pathway to a 
moderate reduction of CH 4 forcing in 2100, but it seems unlikely that the reduction could be 
larger than of the order of 0. 1 W m 2 . 

There is less leverage with N 2 0, whose growth is exceeding all scenarios (Fig. 14b). Major 
780 quantitative gaps remain in our understanding of the nitrogen cycle (Kroeze and Bouwman 

2011) , but fertilizers are clearly a principal cause of N 2 0 growth (Rockmann and Levin 
2005; Park et al 2012). More efficient use of fertilizers could reduce N 2 0 emissions, but 
considering the scale of global agriculture, and the fact that fixed N is an inherent part of feeding 
people, there will be pressure for continued emissions at least comparable to present emissions. 

785 In contrast, agricultural CH 4 emissions are inadvertent and not core to food production. Given 
the current imbalance [emissions exceeding atmospheric losses by about 30% (Prather et al., 

2012) ] and the long N 2 0 atmospheric lifetime (1 16 ± 9 years; Prather et al 2015) it is nearly 
inevitable that N 2 0 will continue to increase this century, even if emissions growth is checked. 
There can be no expectation of an N 2 0 decline that offsets the need to reduce C0 2 . 

790 The Montreal Protocol has been a success in stifling and even reversing the growth of trace 

gases that can destroy ozone and cause global warming (Prather et al 1996; Newman et al 2009). 
Amendments to this protocol to achieve phasedown of additional gases are important, but mainly 
for the objective of limiting the growth of these trace gas climate forcings rather than with an 
expectation of obtaining a large net reduction of climate forcing by MPTGs + OTGs (Fig. 3). 

795 

11. Discussion 

We conclude that the world has already overshot targets for atmospheric temperature and 
greenhouse gas amount required to maintain a safe long-term environment for humanity and 
assure the well-being of young people and future generations. Earth’s paleoclimate history tells 
800 us that, if we wish to avoid locking in multi-meter sea level rise with loss of functionality of 

most coastal cities (Clark et al 2016), our target should be to keep global temperature close to the 
Holocene range, which requires an absolute reduction in current GHG climate forcing and global 
temperature. 

Thus we infer an urgent need for both (1) rapid phasedown of fossil fuel emissions, and (2) 
805 actions that draw down atmospheric C0 2 and, at minimum, eliminate net growth of non-C0 2 
climate forcings. These tasks are formidable and are not now being pursued effectively. 

Although economic and political analysis is outside the scope of this paper, our conclusion 
that the world has already overshot appropriate targets is sufficiently grim to compel us to point 
out that pathways minimizing climate impacts are feasible and have other benefits. The 
810 underlying policy required to spur rapid reduction of fossil fuel emissions is a transparent 

steadily rising carbon fee that makes fossil fuels include their costs to society (Ackennan and 
Stanton 2012; Hsu 2011; Hansen 2014), which encourages energy conservation (reduced 
consumption), energy efficiency, and technology development of carbon-free energy. A rising 
global carbon fee, which could be achieved by agreement of a few major powers (Hsu 2011), is 
815 the crucial underlying policy needed to spur private investment, innovations and consumer 
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choices, but it does not obviate the need for government energy planning, energy efficiency and 
pollution regulations, and support for energy research and development. 

Governments have shown the ability to achieve high rates of emissions reduction, e.g., 

Peters et al (2013) note that Belgium, France and Sweden achieved emission reductions of 4- 
820 5%/year sustained over 10 or more years in response to the oil crisis of 1973. These rates were 
primarily a result of nuclear power build programs, which historically has been the fastest route 
to carbon- free energy (Fig. 2 of Cao et al 2016). Peters et al also note that a continuous shift to 
natural gas led to sustained reductions of 1-2%/year in the UK in the 1970s and in the 2000s, 
2%/year in Denmark in 1990-2000s, and 1.4%/year in the USA since 2005. None of these 
825 examples were aided by the broad economy-wide effect of a rising carbon fee, although high oil 
prices in the 1970s partially simulated that effect. What is needed to achieve rates presently 
demanded by the climate crisis is a combination of a rising carbon fee along with government 
support of technological advances, which has historically received the smallest share of total 
research budgets in OECD countries. 

830 Our scenarios show that, in addition to CO 2 emission phase-out, there must be large CO 2 

extraction from the air and a net halt of growth of non-C02 GHG climate forcings. Success with 
both CO 2 extraction and non-C02 GHG controls requires a major role for developing countries. 
Ancillary benefits of the agricultural and forestry practices needed to achieve C0 2 drawdown, 
such as improved soil fertility, advanced agricultural practices, forest products, and species 
835 preservation, are of interest to all nations. Developed nations have a recognized obligation to 
assist nations that have done little to cause climate change yet suffer some of the largest climate 
impacts. If economic assistance is made partially dependent on verifiable success in carbon 
drawdown and non-CCF mitigation, this will provide incentives that maximize success in carbon 
storage. Similar considerations apply to incentives for reducing trace gas emissions, and, as we 
840 have discussed, some activities such as soil amendments that enhance weathering might be 
designed to support both CO 2 and other GHG drawdown. 

Considering our conclusion that the world has overshot the appropriate target for global 
temperature, and the difficulty and perhaps implausibility of negative emissions scenarios, we 
would be remiss if we did not point out the potential contribution of demand-side mitigation that 
845 can be achieved by individual actions as well as by government policies. Numerous studies (e.g. 
Hedenhus et al 2014; Popp et al 2010) have shown that reduced ruminant meat and dairy 
products is needed to reduce GHG emissions from agriculture, even if technological 
improvements increase food yields per unit fannland. Such climate-beneficial dietary shifts have 
also been linked to co-benefits that include improved sustainability and public health (Bajzelj et 
850 al 2014; Tilman and Clark 2014). Similarly, Working Group 3 of IPCC (2014) finds “robust 

evidence and high agreement” that demand-side measures in the agriculture and land use sectors, 
especially diet shifts, reduced food waste and changes in wood consumption have substantial 
mitigation potential, but they remain under-researched and poorly quantified. 

If rapid emission reductions are initiated soon, it is still possible that at least a large fraction 
855 of required CO 2 extraction can be achieved via relatively natural agricultural and forestry 

practices with other benefits. On the other hand, if large fossil fuel emissions are allowed to 
continue, the scale and cost of industrial CO 2 extraction, occurring in conjunction with a 
deteriorating climate with growing economic effects, may become unmanageable. Simply put, 
the burden placed on young people and future generations may become too heavy to bear. 

860 Appendix A: Additional figures and tables 
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(a) Global Fossil-Fuel C0 2 Annual Emissions 



(b) Global Fossil-Fuel C0 2 Annual Emissions 



Fig. Al. C0 2 emissions from fossil fuel use and cement manufacture, based on data of Boden et al 
(2016) through 2013, with results extended using BP(2016) energy consumption data, (a) is log scale 
and (b) is linear. Growth rates r in (a) for an n year interval from (l+r) n with end-year amount the mean 
for three years to minimize noise. 


Al. Fossil Fuel CO 2 Emissions 

CO 2 emissions from fossil fuels in 2015, based on preliminary data from BP (2016), were only 
slightly higher than in 2014 (Fig. Al). Such slowdowns are common, and usually reflect the 
global economic situation. Given rising global population and the fact that many nations, 
including the soon-to-be-most-populous India, are still at early stages of development, the 
potential exists for continued growth of emissions. Fundamental changes in energy technology 
will be needed if the world is to rapidly change energy course and phase down fossil fuel 
emissions. 


A2. Temperature Data and Analysis Method 

We use the current Goddard Institute for Space Studies global temperature analysis (GISTEMP), 
which is the analysis method described by Hansen et al. [2010] but with updated input data. The 
analysis combines data from three sources: (1) monthly mean meteorological station data of the 
Global Historical Climatology Network (GHCN) described by Peterson and Vose [1997] and 
Menne et al. [2012], (2) monthly mean data from Antarctic research stations of the Scientific 
Committee on Antarctic Research (SCAR), as reported by the SCAR Reference Antarctic Data 
for Environmental Research project (http://www.antarctica.ac.uk/met/READER), and (3) ocean 
surface temperature measurements from the NOAA Extended Reconstructed Sea Surface 
Temperature (ERSST) [Smith et al., 2008; Huang et al., 2015]. 
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12-Montli Ruiuuug Mean Global Temperature 



2020 



Fig. A2a. Global surface temperature (12-month running mean) relative to 1951-1980 in the GISTEMP 
analysis, comparing the current analysis using NOAA ERSST.v4 for sea surface temperature with results 
890 using the prior ERSST.v3b. 

The present analysis uses GHCN.v3.3.0 [Menne et al., 2012] for land data and ERSST.v4 
for sea surface temperature [Huang et al, 2015]. Update from GHCN.v2 used in our 2010 
analysis to GHCN.v3 had negligible effect on global temperature change over the past century 
895 (see graph on http://www.columbia.edu/~mhsl 19/Temperature/GHCN_V3vsV2/). However, the 
adjustments to SST to produce ERSST.v4 have a noticeable effect, especially in the period 1939- 
1945, as shown by the difference between the two data sets (lower graph in Fig. A2a). This 
change is of interest mainly because it increases the magnitude of an already unusual global 
temperature fluctuation in the 1940s, making the 1939-1945 global temperature maximum even 
900 more pronounced than it was ERSST.v3 data. Thompson et al. [2008] show that two natural 
sources of variability, the El Nino/Southern Oscillation and (possibly related) unusual winter 
Arctic warmth associated with advection over high Northern Hemisphere latitudes, partly 
account for global warmth of 1939-1945, and they suggest that the sharp cooling after 1945 is a 
data flaw, due to a rapid change in the mix of data sources (bucket measurements and engine 
905 room intake measurements) and a bias between these that is not fully accounted for. 
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Surface Temperature Anomaly (°C) [1951-80 Base Period] 
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Fig. A2b. Temperature anomalies in three periods relative to 1951-1980 comparing results obtained using 
ERSST.v3b (left column), ERSST.v4 (center column), and their difference (right column). 

Huang et al. [2015] justify the changes made to obtain version 4 of ERSST, the changes 
including more complete input data in ICOADS Release 2.5, buoy SST bias adjustments not 
present in version 3, updated ship SST bias adjustments using Hadley Nighttime Marine Air 
Temperature version 2 (HadNMat2), and revised low-frequency data filling in data sparse 
regions using nearby observations. ERSST.v4 is surely an improvement in the record during the 
past half century when spatial and temporal data coverages are best. On the other hand, the 
largest changes between v3 and v4 are in 1939-1945, coinciding with World War II and changes 
in the mix of data sources. Several hot spots appear in the Southern Hemisphere ocean during 
WW II in the v4 data, and then disappear after the war (Fig. A2b). These hot spots coincide with 
the locations of large SST changes between v3 and v4 (Fig. A2b), which leads us to suspect that 
the magnitude of the 1940s global warming maximum (Fig. 2) is exaggerated; i.e., it is partly 
spurious. We suggest that this warming spike warrants scrutiny in the next version of the SST 
analysis. However, the important point is that these data adjustments and uncertainties are small 
in comparison with the long-term wanning. Adjustments between ERSST. v3b and ERSST.v4 
increase global warming over the period 1950-2015 by about 0.05°C, which is small compared 
with the ~1 °C global warming during that period. The effect of the adjustments on total global 
warming between the beginning of the 20 th century and 2015 is even smaller (Fig. A2a). 
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(a) Total Solar Irradiance (W/m 2 ) 



(b) Sunspot Numbers 
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Fig. A3. Solar irradiance and sunspot number in the era of satellite data. Left scale is the energy passing 
930 through an area perpendicular to Sun-Earth line. Averaged over Earth’s surface the absorbed solar energy 
is ~240 W/m 2 , so the full amplitude of the measured solar variability is -0.25 W/m 2 . 

A3. Solar Irradiance 

Solar irradiance has been measured from satellites since the late 1970s. Fig. A3 is a composite 
935 of several satellite-measured time series. Data through 28 February 2003 are an update of 

Frohlich and Lean (1998) obtained from Physikalisch Meteorologisches Observatorium Davos, 
World Radiation Center. Subsequent update is from University of Colorado Solar Radiation & 
Climate Experiment (SORCE). Historical total solar irradiance reconstruction is available at 
http://lasp.colorado.edu/home/sorce/data/tsi-data/ Data sets are concatenated by matching the 
940 means over the first 12 months of SORCE data. Monthly sunspot numbers support the 

conclusion that the solar irradiance in the current solar cycle is significantly lower than in the 
three preceding solar cycles. 

The magnitude of the change of solar irradiance from the prior solar cycle to the current 
solar cycle is of the order of -0. 1 W/m 2 , which is not negligible but is small compared with 
945 greenhouse gas climate forcing. On the other hand, the variation of solar irradiance from solar 
minimum to solar maximum is of the order of 0.25 W/m 2 , so the high solar irradiance in 201 1 
-2015 contributes to the increase of Earth’s energy imbalance between 2005-2010 and 2010- 
2015. 

950 
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Global Surface Temperature 



Fig. A4. Simulated global temperature with historical climate forcings to 2000 followed by the 
alternative scenario. Historical climate forcings are discussed in the main text. 

955 

A4. Alternative Scenario 

Simulated global temperature for the climate forcings of the “alternative scenario” discussed in 
Section 6 are shown in Fig. A4. Climate model, with sensitivity 3°C for doubled CO 2 , is the 
same as used for Fig. 13. 

960 


Total Greenhouse Gas Forcing Growth Rate 



Fig. A5. Recent growth rate of total GHG climate forcing; see Fig. 9 for individual gases. 

965 A5. Growth Rate of Total GHG Climate Forcing 

In the past several years the growth rate of climate forcing by GHGs has accelerated sharply, in 
contrast to most scenarios, which presumed that the GHG climate forcing would be declining. 
Ozone is not well-mixed, so its changes are not well-measured and are not fully accounted for in 
Fig. 5. However, the effective CH 4 forcing, which is included, includes about half of the 
970 tropospheric O 3 change. 
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Table 

Al. 

Effective Forcing 

(W/m2) Relative to 1850 except Volcanic 

Aerosols 

Year 

co 2 

ch 4 * 

CFCs 

N 2 0 


o 3 $ 

TA+SA 

Volcano 

Solar 


Net 

1850 

0 , 

. 000 

0 , 

. 000 

0 , 

. 000 

0 . 

. 000 

0 . 

. 000 

0 , 

. 000 

-0 . 

. 083 

0 . 

. 000 


-0 . 

, 083 

1860 

0 , 

. 024 

0 , 

. 012 

0 , 

. 000 

0 . 

. 004 

0 . 

. 004 

-0 . 

. 029 

-0 . 

. 106 

0 . 

. 032 


-0 . 

. 059 

1870 

0 , 

. 048 

0 , 

. 025 

0 , 

. 000 

0 . 

. 007 

0 . 

. 009 

-0 . 

. 058 

-0 , 

. 014 

0 . 

. 048 


0 . 

,065 

1880 

0 , 

. 109 

0 , 

. 039 

0 , 

. 000 

0 . 

. 010 

0 . 

. 014 

-0 , 

.097 

-0 . 

. 026 

-0 . 

. 049 


-0 . 

, 001 

1890 

0 , 

. 179 

0 , 

. 054 

0 , 

. 000 

0 . 

. 013 

0 . 

. 018 

-0 . 

.146 

-0 . 

. 900 

-0 . 

. 070 


-0 . 

. 850 

1900 

0 , 

.204 

0 , 

. 073 

0 , 

. 001 

0 . 

. 016 

0 . 

. 023 

-0 . 

.195 

-0 . 

. 040 

-0 . 

.063 


0 . 

, 018 

1910 

0 , 

.287 

0 , 

. 109 

0 , 

. 002 

0 . 

. 020 

0 . 

. 026 

-0 . 

.250 

-0 . 

. 072 

-0 . 

. 043 


0 . 

, 079 

1920 

0 , 

.348 

0 , 

. 150 

0 , 

. 003 

0 . 

. 027 

0 . 

. 032 

-0 , 

.307 

-0 , 

.215 

-0 . 

. 016 


0 . 

. 022 

1930 

0 , 

. 425 

0 , 

.194 

0 , 

. 004 

0 . 

. 035 

0 . 

. 036 

-0 , 

.364 

-0 , 

.143 

0 . 

. 014 


0 . 

,200 

1940 

0 , 

.494 

0 , 

.232 

0 , 

. 005 

0 . 

. 041 

0 . 

. 045 

-0 , 

. 424 

-0 , 

. 073 

0 . 

. 037 


0 . 

,356 

1950 

0 , 

.495 

0 , 

.274 

0 , 

. 009 

0 . 

. 049 

0 . 

. 056 

-0 . 

.484 

-0 . 

.066 

0 . 

. 055 


0 . 

.387 

1960 

0 , 

.599 

0 , 

.342 

0 , 

. 027 

0 . 

. 057 

0 . 

. 078 

-0 . 

. 621 

-0 . 

. 106 

0 . 

. 102 


0 . 

, 478 

1970 

0 , 

.748 

0 , 

. 433 

0 , 

. 076 

0 . 

. 071 

0 . 

.097 

-0 , 

.742 

-0 , 

.381 

0 . 

.093 


0 . 

.395 

1980 

0 , 

. 976 

0 , 

.532 

0 , 

. 185 

0 . 

.091 

0 . 

.115 

-0 . 

. 907 

-0 , 

. 108 

0 . 

.169 


1 . 

. 054 

1990 

1 . 

.227 

0 , 

. 618 

0 , 

.303 

0 . 

.118 

0 . 

.117 

-0 . 

. 997 

-0 , 

.141 

0 . 

. 154 


1 . 

,399 

2000 

1 . 

.464 

0 , 

. 651 

0 , 

.347 

0 . 

.141 

0 . 

.117 

-1 . 

. 084 

-0 . 

. 048 

0 . 

. 173 


1 . 

.761 

2005 

1 . 

. 619 

0 , 

. 651 

0 , 

.356 

0 . 

. 153 

0 . 

. 123 

-1 . 

. 125 

-0 . 

. 079 

0 . 

. 019 


1 . 

,716 

2010 

1 . 

.766 

0 , 

. 665 

0 , 

.364 

0 . 

.167 

0 . 

. 129 

-1 . 

.163 

-0 . 

. 082 

0 . 

. 028 


1 . 

, 874 

2015 

1 . 

. 927 

0 , 

. 684 

0 , 

.373 

0 . 

. 183 

0 . 

. 129 

-1 , 

.199 

-0 , 

. 100 

0 . 

. 137 


2 . 

. 134 

# CH 4 : CH 4 -induced changes of tropospheric 0 3 and stratospheric H 2 0 are included. 
$ 0 3 half of tropospheric 0 3 forcing + stratospheric 0 3 forcing from IPCC (2013) 
Annual data are available in a longer version of the table available at 

http://www.columbia.edu/~mhsll9/Burden/ . 


Table 

A2 . Effective 

Forcing 

(W/m 2 ) Relative 

to 1850 

except Volcanic 

Aerosols 

Year 

Q 

O 

# ch 4 

CFCs 

N 2 0 

$ o 3 

TA+SA 

Volcano 

Solar 

Net 

2016 

1 . 942 

0 . 654 

0.367 

0 .180 

0 . 130 

-1.207 

-0 . 100 

0 .097 

2 .062 

2020 

2 . 074 

0 . 658 

0.373 

0 .189 

0 . 130 

-1.234 

-0 . 100 

-0 . 008 

2 . 082 

2030 

2.347 

0 . 663 

0.343 

0.212 

0 . 130 

-1.296 

-1 . 057 

-0 . 008 

1.335 

2040 

2.580 

0 . 688 

0.301 

0.238 

0 . 123 

-1.350 

-0 . 100 

0 . 027 

2.507 

2050 

2 . 803 

0.717 

0.267 

0.271 

0 .117 

-1.396 

-0 . 100 

0 .062 

2.741 

2060 

3 . 017 

0.740 

0.243 

0.302 

0 .111 

-1 . 433 

-1.208 

0 .097 

1 . 870 

2070 

3.222 

0.753 

0.229 

0.337 

0 . 105 

-1.462 

-0 . 100 

0 . 132 

3.215 

2080 

3 . 421 

0.741 

0.215 

0.367 

0 .098 

-1 .484 

-0 . 100 

0 .167 

3 . 425 

2090 

3 . 614 

0 . 676 

0 .199 

0 .401 

0 .091 

-1.495 

-1.240 

0 .167 

2 .413 

2100 

3 . 801 

0.580 

0 .191 

0 . 428 

0 . 085 

-1.500 

-0 . 100 

0 .167 

3 . 652 


# CH 4 : CH 4 -induced changes of tropospheric 0 3 and stratospheric H 2 0 are included 
? 0 3 : Half of tropospheric 0 3 forcing + stratospheric 0 3 forcing from IPCC 2013 
Annual data are available in a longer version of the table available at 

http://www.columbia.edu/~mhsll9/Burden/ . 
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